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1 . 


INTRODUCTION 


Manmade satellites have provided data for some of the more 
important scientific contributions of our century, and quite justifi- 
bly mankind should have great admiration for those whose ingenuity and 
imagination have made them possible. Yet long before the Wright 
Brothers accomplished their historic flight at Kitty Hawk, long before 
Fulton engineered the steamboat, and even before James Watt unraveled 
the principles behind the steam engine, a man named Isaac Newton 
considered the possibilities of satellites and outlined the conditions 
necessary to make them possible. Three hundred years ago, in 
Principia Mathematica , Newton discussed the conditions under which a 
projectile must be launched to follow a specific orbit and laid the 
precise mathematics of orbital motion. Newton foresaw rocket propul- 
sion and was able to state the theory behind it. Over three centuries 
ago, Isaac Newton understood the mathematic principles that are 
readily accepted and applied in today's world to propel spacecraft 
into the heavens. 

Since the beginning of the space flight program of the United 
States in 1957, many various unmanned spacecraft have been launched 
successfully and have ultimately revealed many new horizons that were 
previously only imagined. With the launch of TIROS 1, the first full- 
time meteorological satellite, on April 1, 1960, it became quite 
evident that the satellite and the data generated from it proved 

The citations on this and following pages follow the style of 
the Journal of Applied Meteorology . 


worthwhile and exciting when applied to meteorological problems. The 
TIROS series extended man's understanding of the formation and dis- 
tribution of various cloud systems and produced data on global cloud 
cover which provided valuable information on global heat balance and 
water vapor distribution. Furthermore, the TIROS series provided the 
meteorologist with almost real time photographs of weather systems 
and permitted early detection of hurricanes, thereby permitting more 
accurate forecasts. 

As technology rapidly progressed and the meteorological satellite 
program grew, vast numbers of special applications were investigated. 
Rainfall estimation was one of these. 
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2 . 


BACKGROUND 


Lethbridge (1967), at Pennsylvania State University, conducted 
one of the earliest studies comparing local rainfall with image 
brightness and infrared (IR) temperatures. As technology advanced 
over the next 6 years and systems for data transmission and acquisi- 
tion improved, more detailed studies were begun. Pioneers in these 
efforts were Simpson and Woodley (1971), Sikdar (1972), Martin and 
Scherer (1973), and Follansbee (1973). 

Martin and Scherer (1973) furnish a fairly complete review of 
methods to estimate rainfall (up to 1973) from either visible (VIS) or 
IR satellite imagery, Follansbee (1973) primarily deals with rain- 
producing clouds in his empirical approach to areal daily average rain- 
fall estimations. 

More recent applications of rainfall estimation have been con- 
ducted by Allison ^ al. (1974) using NIMBUS 5 electronically scanning 
microwave radiometer data over oceans; and by Griffith ^ (1976) 

using Applications Technology Satellite (ATS-3) imagery over southern 
Florida. 

The first rainfall estimation studies that used geosynchronous 
satellite imagery were conducted by Follansbee and Oliver (1975) and 
Oliver and Scofield (1976). As a result of those studies, a scheme for 
point approximation of convective rainfall was developed by Scofield 
and Oliver (1977) using Geostationary Operational Environmental Satel- 
lite (GOES) VIS and enhanced infrared (EIR) imagery. Recently, Scofield 
(1978a, 1978b) emphasized the practical applications of the GOES-1 VIS 
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and EIR imagery to detect and estimate rainfall from thunderstorm- 
produced flash floods by expanding his point approximation scheme to 
include areal coverage. Since GOES VIS and IR imagery are provided to 
the user every 30 min during the system’s normal operation and can be 
adjusted to provide data at 5-min intervals during rapid scan opera- 
tions (Corbell ^ 1976), a careful severe weather (e.g., flash 

flood) watch can be maintained over vast areal expanses. 

The May 1978 issue of the Bulletin of the American Meteorological 
Society presented a statement of concern by the American Meteorological 
Society (1978) regarding actions that must be taken to prevent loss of 
life in flash floods. Their position makes it clear that much more 
must be learned about rapid convective rainfall if we are to provide 
an objective input into the decision making process regarding flash 
flood forecasting and the associated warnings to the public. Mogil ^ 
al . (1978) outlined the current flash flood warning program of the 
National Weather Service and stressed that additional research is 
necessary in various data-gathering operations, one of these being the 
development of techniques for rainfall estimation based on satellite 
data . 

At present, the most accurate rainfall measurements are those 
obtained from a close network of raingages. It is obvious that rain- 
gages cannot be placed everywhere, let alone in a sufficiently dense 
network that would be economically feasible. Rugged mountainous ter- 
rain, remote interior areas on continents, and vast ocean areas are 
other apparent factors limiting such a raingage network. Even over 
flat terrain, a close network of raingages may not be representative 
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if the rainfall is from convective activity. Therefore, alternative 
methods for approximating rainfall must be sought. 

The meteorological radar system is well established and receives 
much attention from forecasters and hydrologists for precipitation 
estimates. Though it is less accurate than a close network of rain- 
gages, the radar system is economically feasible, in most cases, and 
gives a very good quantitative estimate over large areas. However, 
the weather radar is chiefly limited to land based operations for 
rainfall estimates. Furthermore, the employment and maintenance of a 
complex radar system is dependent upon the user's technology and 
operating capital, thereby limiting its use to the developed nations. 
Therefore, with the limitations placed on both the raingage networks 
and weather radars, satellite imagery for rainfall prediction appears 
to be a viable solution to the problem. 

Reynolds and Vender Haar (1973) compared the reflected solar 
radiation from ATS-3 imagery with cloud heights determined by using 
radar. Their study showed a correlation of 0.88, which they feel is 
significant enough to support the argument that the cumulus height and 
growth may be inferred from geostationary satellites. Studies by 
Woodley et ad. (1972), Griffith and Woodley (1973) and Griffith ^t al . 
(1976) substantiate the same premise that the brighter the cloud on 
visible imagery, the higher the cloud and the heavier the rainfall. 

This principle-- the higher (colder) the clouds, the heavier the 
rainfall--was the very basis for later work. 

Following the GARP (Global Atmospheric Research Program) Atlantic 
Tropical Experiment (GATE), Griffith ^ (1976) reported on their 
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findings from ATS-3 imagery to formulate a procedure for approximating 
rainfall during daylight hours. The ATS-3 data were digitized by use 
of a densitometer and later calibrated via gage-adjusted 10-cm radar. 
They suggested that higher resolution data (e.g., digital GOES data) 
would have improved the accuracy of their results by 20 to 307o. In 
a recent paper, Griffith ^ (1978) discussed the feasibility of 

real time flash flood prediction based on the methodology from their 
GATE investigations. They presented cumulative rainfall results from 
the Big Thompson storm (July 31 and August 1, 1976) from GOES IR 

imagery. Although their approach was somewhat imprecise, they feel 
that such data could be used in near real time. Once again this 
demonstrates that digital data are superior to the imagery and should 
be used when available. 

Restraints imposed on a scheme for rainfall estimation from 
digital GOES data by use of an objective technique would be fairly 
rigorous. First, the satellite data must be "navigated" (rectified) 
and then reduced to a suitable spatial resolution. Next, the cloud 
system of the individual cloud must be isolated and tracked so that 
rainfall amounts can be computed per grid square over a time interval. 
Some of the latest work has been accomplished by Scofield (1978b), who 
demonstrates that his methods using GOES VIS and IR imagery can be 
applied to flash flood problems. The scheme is relatively simple and 
is not burdened by the problems of navigation and spatial reduction. 

His method does, however, rely heavily upon the experience level of the 
individual interpreting the satellite imagery. 

The most important, obvious step when handling satellite imagery 


6 



is to accurately grid (navigate) the data to provide maximinn resolu- 
tion. Roth (1969) decribes an automated technique for the rectifica- 
tion and mapping of meteorological satellite imagery. His method is 
subjective, biased toward data selection; whereas. Wolf et al, (1972) 
and Cherne (1974) improved on the technique Roth developed and 
navigated the satellite Imagery while devising a technique to deter- 
mine cloud motions. They basically relate the location of the satel- 
lite in time and space to the actual surface of the earth and then 
compare the actual surface and the map projection. Cherne (1974) 
also describes the immense volumes of data that severely tax even 
automated methods of data processing. 

Hawkins (1977) agrees with Cherne and stresses that the large 
amounts of data place limitations on all aspects of its use. There- 
fore, data amounts must be reduced as efficiently as possible while 
retaining maximum resolution. Another way to minimize the vast 
amounts of data and retain the maximum resolution would be to observe 
a significantly smaller areal expanse. The GOES imagery system does 
this by allowing the user to observe a smaller areal expanse without 
sacrificing resolution. 

Isolating and tracking areas of interest could simply be 
observed from the imagery in a subjective manner (like Scofield and 

Oliver, 1977) or by an automated, objective process, whereby exceeding 
a specified threshold value would trigger a signal for a closer 
weather watch in that sector. Since an experienced analyst in satel- 
lite photo interpretation can take a glimpse at a few pictures of some 
area of the globe and make some s umm arizing statements packed with 
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meteorological information, it seems impossible that an objective, 
automatic technique could handle all the special cases (Hawkins, 1977). 
Though it seems that a compromise between the objective and subjective 
analysis methods will be necessary for some time, the potential value 
of an objective technique for rainfall estimation is great. 

In light of the above, it is hoped that this investigation will 
be another step toward the formulation of an objective technique for 
rainfall estimation and will become one that v/ould be used as a 
supplemental tool by meteorologists and hydrologists responsible for 
flash flood prediction. 
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3. 


THE GOES-1 SYSTEM 


The main goal of the GOES satellite program is to increase 
knowledge and understanding of earth's atmosphere and the processes 
that take place within it by observing the life cycle and motion of 
s 3 moptic systems, tropical cyclones, severe thunderstorms, and other 
related atmospheric phenomena. Another objective of the program is to 
develop a national and international environmental network that can 
obtain, process, and transmit the GOES data and weather related 
warnings in real time (McKowan, 1977). 

GOES-1 was launched October 16, 1975 and arrived on station 
December 18, 1975, Its position is 35,800 km above the equator at 
75“W. It travels at a speed of about 11,000 km h ^ in the same 
direction as the earth's rotation on its axis. At this speed and 
altitude the satellite remains relatively fixed at its position in 
space with respect to the surface of earth, and is thus termed a 
geosynchronous satellite. The satellite maintains a spin rate of 100 
rpm as its spin axis remains parallel to the earth's polar axis. 

This enables the spacecraft to maintain a perpendicular orientation to 
the earth equatorial plane (Corbell ^ 1976). 

The visible infrared spin-scan radiometer (VISSR) is capable of 
day and night mapping of the earth and the clouds in earth's atmo- 
sphere because of the satellite's synchronous spin- stablized geo- 
stationary orbit. The VISSR scans from west to east (left to right) 
in eight identical VIS channels (0. 55-0 . 70 yum band) and two redundant 
IR channels (10.5-12.6 yum band). The VIS and IR data sets are 
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produced simultaneously (see Fig. 1) and a full disc of earth as 
viewed by the IR sensor would consist of 1,821 latitude steps (scan 
lines) with 3,822 points per step. Thus, the IR^ sensor detects 
6,959,862 separate IR picture elements (pixels). The corresponding 
resolution is 4 n mi by 2 n mi at satellite subpoint. Furthermore, 
the associated mode A VIS image, taken simultaneously with the IR 
image, would consist of 14,568 scan lines and have 15,288 points per 
line. For a full disc of earth that would amount to more than 222 
million pixels with a subpoint resolution of % n mi by % n mi 
(McKowan, 1977). 

At this point, it is noteworthy to pause and consider the 
enormous amounts of data that are generated by the GOES system, as 
well as the methods of data retrieval and redistribution. Since a 
full disc IR image amounts to almost 7 million pixels, it seems fairly 
obvious that some data reduction and sectorizing must take place to 
allow for near real time transmission to the Satellite Field Services 
Stations (SFSS) and the other users. The unprocessed VISSR data are 
transmitted to the Command and Data Acquisition (CDA) station’s 18.3 m 
antenna, at Wallops Island, Virginia. The CDA station processes the 
raw data through a synchronizer/data buffer which separates the VIS 
and IR data. The VIS data are then transmitted back to the satellite, 
while the lower resolution IR data are processed through computers 
for analog transmission via 3 kHz telephone lines direct to the 
Central Data Distribution Facilities (CDDF) at Marlow Heights, Mary- 
land, The VIS data, received from the satellite, are then relayed 
to the CDDF via a microwave link from the National Environmental 
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VISIBLE SENSOR LINES 



:ture data showing the relationship between infrared and 


Satellite Service (NESS) at Suit land, Maryland, The data are then 
formated and sectorized by CDDF to allow telephone line transmission 
to the SFSS's. Fig. 2 illustrates the data flow and depicts the 
transmission from location to location and the arrival at the SFSS's 
and the other users. Table 1 represents the locations of SFSS's 
and their co-located National Weather Service (NWS) units. 

Table 1. Current SFSS's in operation (from Corbell ^ al., 1976). 




SFSS 

Co-located NWS Units 


Washington, D.C. Washington, D.C. Weather Service Forecast 

Office (WSFO) 


Miami, Florida National Hurricane Center, Miami 

Kansas City, Missouri National Severe Storms Forecast Center 

Kansas City, Missouri 

San Francisco, California San Francisco WSFO and Eastern Pacific 

Hurricane Center 

Honolulu, Hawaii Honolulu WSFO and Central Pacific 

Hurricane Center 

Anchorage, Alaska Anchorage WSFO, Alaska 


The Wallops Island CDA station converts the radiance measurements 
to equivalent blackbody temperatures (the calibration process). The 
CDDF sectorizes the digital data and transmits them to the users. The 
256 digital count values, assigned during the calibration process, 
represent specific ranges of temperature. Each count represents a 0.5K 
change from 330K to 242K and a l.OK change from 242K to 163K (Fig. 3). 
Since it is difficult to recognize significant features in IR imagery 
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Fig. 2. GOES data flow diagram (after Corbell et al., 1976). 


Fig. 3, Standard IR calibration of digital count values to their 
associated equivalent blackbody temperatures (from Corbell 1976). 
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when temperature gradients are small, enhancement’s main purpose is to 
provide contrast between fog and haze, thunderstorm tops, ocean 
current boundaries, and terrain features. The enhancement curves that 
are utilized with IR imagery furnish field forecasters with another 
tool by which an improved weather forecast, watch, and/or warning can 
be disseminated to the public. The assorted enhancement curves avail- 
able to SFSS's are outlined by Corbell ^ (1976). 

The enhancement curve that receives the most attention for 
convective precipitation estimation is the Mg curve because it allows 
better definition of very cold domes (Fig, 4). One of the drawbacks, 
however, is that specific temperatures below -65° C cannot be obtained 
from the IR imagery. It does, however, progress from black to white 
gradually rather than producing a complete white-out at temperatures 
colder than -65° C. This provides the user with a good idea just where 
the coldest tops are and consequently locates the areas of maximum 
precipitation. 

Scofield and Oliver (1977) have used the Mg curve as a basis for 
their rainfall estimation studies. This investigation uses the con- 
cept of enhancement as the foundation for this investigation. 

Inasmuch as this research primarily deals with digital GOES 
data, further discussion will be limited to details related to this 
investigation. Information concerning subsystems of the spacecraft, 
archival of the data, and other details can be found in Corbell e^ al . 
(1976) or McKowan (1977). 
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COUNT VALUE (INPUT) 


WHITE 


BLACK 
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Number 


— f 000 
-165 °F 
-103 °C 
164 °K 


TEMPERATURE 


° C Temperature 
to 

4*58.8 to +28.8 
+28.8 to + 6.8 
+ 6.8 to “31.2 
-32.2 to -41.2 
-42.2 to -52.2 
-53.2 to -58.2 
-59.2 to -62.2 
-63.2 to -80.2 


-80.2 to -109.0 


Comments 

Reason for enhancement 

Little or no useful data 
Low level/ sea surface 

Middle level no enhancement 

First level contour (Med Grey) 
Tstm enhancement (Light Grey) 
Tstm enhancement (Dark Grey) 
Tstm enhancement (Black) 
Overshooting tops enhancement 
(gradual increase from Black 
to White) 

(White) 


Fig. 4. Mg enhancement curve (after Corbell et al., 1976). 
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4. 


DATA 


In order to examine the feasibility of a convective rainfall 

V 

approximation scheme from digital GOES-1 data, two primary data bases 
are required, as well as some supplemental, supportive data. First, 
it is necessary to construct a ground truth for comparisons with the 
satellite rainfall estimations. For this part of the study, data from 
the National Climatic Center (NCC) publications for May 1978 were used 
in order to determine raingage precipitation amounts for the 24-h and 
hourly time segments. Second, and most important, the digital GOES-1 
data were obtained from the Man Computer Interactive Data Access 
System (McIDAS) at the University of Wisconsin. McIDAS receives the 
satellite data directly from the raw GOES data base. 

a. Raingage Data 

The NCC data were carefully extracted from the May 1978 publica- 
tions, Climatological Data (Vol. 83, No. 5) contains the 24-h precip- 
itation amounts for stations in the 480 km by 480 km square centered 
on College Station (CLL), Texas. The Hourly Precipitation Data (Vol. 
28, No. 5) contains rainfall rates for the hourly time intervals of 
interest. National Weather Service radar summaries and Texas A&M 
University (TAMU) digital radar information supplemented the analysis 
of the raingage amounts by enabling a fairly accurate positioning of 
precipitation areas versus precipitation-free areas. The two forms of 
radar data indicated speed and direction of movement of the system and 
individual cells, as well as a clarification of rain intensity and 
location . 
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b 


Digital GOES-1 Data 


As mentioned earlier, a single satellite image consists of 
millions of bytes of information, and the data must be handled care- 
fully to reduce them to specific areas of interest while maintaining 
maximum resolution. The McIDAS data used in this investigation were 
produced with VIS mode data of 1 km by 1 km resolution and IR mode 
data with 4 km by 4 km resolution. The VIS data were not used because 
the night terminator moved across the study area during the first hour 
of investigation and rendered the data of little value for rainfall 
estimations due to lack of visible light. An investigation of VIS 
data during the daytime hours would be valuable because the VIS data 
have superior resolution that would allow even better raingage and 
radar comparisons. 

The IR data, which provide nearly continous viewing of earth, 
have 4 km resolution by virtue of the method with which the McIDAS 
handles the data. Whereas, the GOES system records an IR pixel in a 
4 km (horizontally) by 8 km (vertically) resolution (McKowan, 1977), 
the McIDAS simply divides the pixel into two 4 km by 4 km pixels and 
displays them over two adjacent sensor lines (Young, 1978). The 
resolution, therefore, is implied as 4 km, but in reality it is not 
because the VISSR integrates the signal over the 8 kra interval. 
However, for data handling, the 4 km pixel appears to be in a less 
confusing and more usable form. 

The GOES data were extracted from the nine-track McIDAS tapes by 
using a computer programming technique developed by Henderson (1979) . 
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c. 


Supplemental Data 


Supportive data for this investigation were gathered from 
several different areas. The TAMU digital weather radar provided data 
for the same time frame as this study. The VIL (vertical integration 
of liquid water content) maps, produced by McAnelly (1979), were used 
to fine tune the raingage analysis, and then later were compared to 
the digital GOES-1 data. The rawinsonde and some synoptic data that 
were used in this investigation came from NASA's AVE VII (Atmospheric 
Variability Experiment) experiment^ (Davis et al. y 1978). The remain- 
der of the NWS synoptic data, as well as GOES-1 images. Weather Wire, 
and Service "C" teletype data were gathered from the TAMU Department 
of Meteorology archive. 


1. Atmospheric Sciences Division, Space Sciences Laboratory, NASA 
Marshall Space Flight Center, Alabama 35812 
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5. 


RESULTS 


This investigation was conducted over the region of Texas 
depicted by Fig. 5. College Station (TAMU) is positioned in the 
center of the two concentric boxes. The large box is 480 fcm by 480 
km, and the smaller box is 240 km by 240 km. The raingage analyses 
were performed over the larger box to ensure continuity, but will be 
presented here only for the smaller box to allow maximum resolution 
for the comparisons with the digital GOES-1 data. The GOES data were 
analyzed for rainfall within the 240 km by 240 km box. 

a. The Synoptic Situation 

At 1820 CDT (2320 GMT) on Tuesday, May 2, 1978, the National 
Severe Storm Forecast Center at Kansas City issued Tornado Watch 
Number 97 for portions of central and southeastern Texas and adjacent 
coastal waters. The watch, valid from 1900 CDT until 0100 CDT (0000-' 
0600 GMT), was issued because a strong upper level low was tracking 
across the south central United States and triggering severe thunder- 
storms and associated phenomena as it moved slowly eastward above a 
somewhat stationary frontal system (Fig. 6a), The warm front 
stretched across Texas from Del Rio to San Antonio to Houston and on 
eastward toward Louisiana, This sluggishly moving frontal system and 
the excess moisture that was stacked from low to mid levels in the 
troposphere, would continue to interact with the upper level flow 
(Fig. 6b), and thus maintain the potentially severe weather pattern. 

As the upper level low moved eastward, numerous warnings of 


21 




Fig. 5, Area under investigation (small box is 240 km 
by 240 km and the larger box is 480 km by 480 km) with CLL 
located at the center. 
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Fig, 6b, 300 mb upper level wind at 0000 GMT, May 3, 1978. 
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severe thunderstorms were issued by regional NWS Offices (NWSO) via 
the weather wire. A severe weather statement issued by the NWSO at 
Houston at 2005 CDT (0105 GMT) noted that the tornado watch was still 
in effect for most of southeast Texas and that the thunderstorm 
activity nearest to East Texas was just west of Somerville Dam. The 
heaviest thunderstorm activity occurring at 2035 CDT (0135 GMT) was in 
the vicinity of the line echo wave pattern (LEWP) located in South 
Central Texas, as illustrated by Fig. 7. The NWS reported rains in 
excess of 50 mm h ^ from those storms in and near the LEWP and further 
reported that all precipitation was moving to the east at 13 m s 

At 2115 CDT (0215 GMT) the NWSO at San Antonio issued a flash 
flood watch for the remainder of the night for the area along and east 
of a line from Corpus Christi to Victoria to College Station. The 
forecasters believed that the heavy thunderstorm line (LEWP) moving 
into that region would produce rains in excess of 50 mm with the 
potential for some localized rainfall amounting to near 125 mm. By 
2135 CDT (0235 GMT), the NWS reported a line of heavy thunderstorms, 

24 km wide, just west of College Station (Fig, 8). This line extended 
southward into South Central Texas and was reportedly moving to the 
southeast at 13 m s Other thundershowers and isolated heavy 

thunderstorms were observed 40 km northeast of College Station and 
also in the area of Lakes Conroe and Livingston. At 2230 CDT (0330 
GMT), however, the NWSO at Austin cancelled the tornado and flash 
flood watches for much of central Texas because the heaviest activity 
had moved into East Texas where the tornado and flash flood watches 
still remained in effect. 
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As the evening progressed the thunderstorm activity pressed 
further eastward and began to dissipate rapidly while leaving only 
scattered thunderstorms across eastern Texas and Jthe upper coastal 
areas. At 0145 CDT (0645 GMT), Wednesday, May 3, 1978, the NWS 
cancelled all the watches in Texas because the system had decreased 
sufficiently in both areal coverage and intensity, and no longer posed 
a threat for heavy rain or flash flooding. 

b. The Raingage Analysis 

Careful analyses were performed on the reported raingage amounts 
collected on the evening of May 2, 1978, across central and eastern 
Texas. These are presented over the next few pages. Fig. 9 repre- 
sents the 24-h (02/1200 GMT - 03/1200 GMT) rainfall analysis for the 
area of investigation within the 240 km by 240 km box centered on 
TAMU (CLL) . Although this portion of the investigation closely 
examines a 3-h segment (0100-0400 GMT) within the 24-h period, it is 
noteworthy to point out that much of the precipitation fell within a 
3- to 5-h period, and in one instance within 1-h. Table 2 depicts the 
amounts of rainfall reported and the associated time intervals within 
which precipitation occurred. 

To ensure continuity with the 24-h data, each hourly segment of 
precipitation was analyzed over the 480 km by 480 km box and then 
presented here on the smaller, 240 km by 240 km box. A technique 
using the NWS radar summaries, TAMU's digital radar data, NCC s 
hourly precipitation data, graphical addition and subtraction of 
isohyet fields of rainfall, and the author's 15 years of observing and 
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Fig. 9, 24-h rainfall analysis (02/1200 GMT - 03/1200 GMT) in 

nnn for a 240 km by 240 km box centered on CLL. 
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Table 2. Hourly rainfall amounts and 3-h/24-h total rainfall in mm 


Location 

First 

hour 

Second 

hour 

Third 

hour 

3-h total 

24-h total 

Belton Dam 

4 

7 

0 

11 

38 

Conroe 

0 

0 

0 

0 

10 

Eagle Lake 

0 

0 

5 

5 

10 

Groesbeck 

3 

36 

5 

44 

48 

Houston-WSMO 

0 

3 

0 

3 

10 

Houston-Addicks 

1 

1 

13 

15 

21 

Houston- Sat suma 

0 

1 

0 

1 

20 

Jewett 

0 

8 

16 

24 

28 

Lexington 

3 

33 

5 

41 

48 

Love lady 

0 

0 

14 

14 

21 

Post Oak School 

2 

1 

17 

20 

26 

Red Rock 

0 

20 

0 

20 

25 

Richmond 

0 

1 

4 

5 

16 

Somerville Dam 

3 

0 

13 

16 

20 

Stillhouse Hollow Dam 

4 

4 

0 

8 

25 

Waco Dam 

12 

4 

1 

17 

25 

Washington State Park 

0 

3 

8 

11 

18 

Wheelock 

4 

11 

3 

18 

20 
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forecasting experience was used to analyze Figs. 10-13. Figs. 10, 11, 
and 12 represent each hour of rainfall, whereas, Fig. 13 is the total 
rainfall for the 3-h interval. 

c. The Satellite Analysis 

The digital GOES-1 data were recovered from nine-track McIDAS 
tapes by using a computer programming technique by Henderson (1979) 
and subsequently made available for this rainfall estimation study. 

The navigated, digital count values are plotted on a 76 cm by 76 cm 
grid that is representative of the 240 km by 240 km box centered on 
TAMU. The satellite data obtained from McIDAS contained information 
for each 15 min interval from 03/0015 GMT until 03/0345 GMT for a 
total of 15 digital count value maps (DCVM) . The first two maps were 
omitted since little activity of significance was occurring. The 
remaining 13 DCVMs were used in this study and are presented for 
comparison with the images based on the Mg enhancement curve (Mg EIR) 
for 0100, 0200, 0230, 0300, and 0400 GMT over the next few pages. 

Naturally, digital counts ran from low (63) to high (222) count 
values over the computer gridded output maps. However, to make 
analysis less confusing and to enable simplified presentation of the 
data, only selected isopleths of the digital counts are shown. 
Furthermore, since the digital count values provide little meaning as 
to temperature without a conversion method. Fig. 3 (p. 15) was pro- 
vided to relate the digital count value to its equivalent blackbody 
temperature. 

The first step before beginning any satellite-derived rainfall 
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Fig. 10. First hour (0100-0200 GMT) rainfall analysis in mm. 
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Fig, 11. Second hour (0200-0300 GMT) rainfall analysis in mm. 
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Fig. 12. Third hour (0300-0400 GMT) rainfall analysis in mm. 
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Fig. 13. 3-h (0100-0400 GMT) total rainfall in mm. 
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approximations was to draw isopleths of the digital count values. 

After that, a digital count/ observed rainfall table was constructed to 
better understand the variations of digital count values to rainfall 
amounts. To remain as objective as possible during this stage of the 
research, only rainfall reporting stations within the study area and 
under convective activity were used. Table 3 illustrates a portion of 
the original table prepared for the second hour in the northwest 
quadrant of study. The given quarter-hour digital count for each 

Table 3. Digital count /observed rainfall^for the northwest quadrant 
during the second hour (rainfall in mm h ). 


Time (GMT) 

Waco Dam 

Groesbeck 

Belton Dam 

Stillhouse Wheelock 
Hollow Dam 

0200 

215 

212 

211 

211 

209 

0215 

211 

215 

211 

211 

211 

0230 

217 

219 

218 

217 

219 

0245 

209 

217 

209 

210 

214 

0300 

206 

214 

208 

199 

215 

Rainfall 

4 

36 

7 

4 

11 

station in 

question was 

recorded. 

Then the variation of digital 

counts 

was empirically compared 

to actual 

precipitation 

amounts observed 

. By 


experimenting with the digital counts and observed rainfall amounts, 
various count values were tested to determine the best thresholds for 
rainfall estimates. The digital count 208 proved to give the best 
estimates when compared with observed cases. Therefore, this study 
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is based primarily on 208 and higher digital count values. 

Of particular interest is the fact that the 208 count value is 
the overshooting tops enhancement (segment 8) of the Mg enhancement 
curve shown earlier in Fig. 4 (p. 17). That portion of the % curve 
gradually increases from black (count 208) to white (count 225) and 
only the very coldest tops are easily seen on the imagery. This will 
be illustrated further when the 0230 GMT imagery and the DCVM are 
compared. Furthermore, the availablity of digital count values rather 
than an image threshold appears to indicate that they are important to 
consider and may very well be the most accurate means to adjust, and 
thus improve the operational forecast for rainfall amounts in select 
topographic regions. The conclusions section will discuss this 
concept further. 

Fig. 14 for the 3rd of May at 0045 GMT, depicts the first 
overshooting tops (count 208-225) moving into the area of study. The 
0100 GMT GOES-1 Mg EIR image was reproduced and is presented in Fig. 

15. The extent of the storm and the fact that the high, cold cloud 
tops are widespread should be noted. As the system pushed further 
eastward (Figs. 16a-d) and the line of thunderstorms moved into the 
western part of the area, the first of the highest digital counts 
(220-222) was observed (Fig. 16d) . Fig. 17, the 0200 GMT satellite 
image, shows the area where the highest digital counts (coldest 
temperatures) were observed. They do not, however, provide informa- 
tion that would be usable for high resolution precipitation evalua- 
tions. The DCVMs for the next hour (Figs. 18a&b; 20a&b) are the most 
interesting and revealing in this investigation- The 0200 GMT and 
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Fig. 14. Digital count value map (DCVM) at 03/0045 GMT. 
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Fig. 16a. DCVM at 03/0100 GMT. 
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Fig. 16b. DCVM at 03/0115 GMT. 
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Fig. 16c. DCVM at 03/0130 GMT. 
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Fig. 16d. DCVM at 03/0145 GMT. 
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Fig. 18a. DCVM at 03/0200 GMT. 
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Fig. 18b. DCVM at 03/0215 GMT. 
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Fig. 19. 0230 GMT GOES-1 Mg EIR image. 
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Fig. 20a. DCVM at 03/0230 GMT. 
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Fig. 20b. DCVM at 03/0245 GMT. 
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0215 GMT DCVMs depict the advance of the line of thunderstorms and 
their gradual strengthening prior to 0230 GMT. The 0230 GMT GOES 
image (Fig. 19) and DCVM (Fig. 20a) represent the period of maximum 
horizontal and vertical growth during the 3-h time frame. The second 
hour was the period of maximum rainfall also. 

Fig. 19 best illustrates the concept introduced earlier, that 
only the very coldest tops are easily observed on the imagery. Note 
the area of bright white pixels in and on either side of the maximum 
activity in the thunderstorm line. A quick examination of Fig. 20a 
shows the region of high digital count values that are otherwise 
camouflaged in segment 8 (count 208-225) of the Mg enhancement curve. 
The 0245 GMT DCVM (Fig. 20b) marks the collapse of the strongest 
period that was occurring in the previous DCVM and denotes the initial 
stages of decay of the system as a whole. Fig. 20b corresponds to the 
end of the hour of heaviest precipitation (compare with Fig. 11, p. 

32) . 

The 0300 GMT satellite image (Fig. 21) substantiates the thesis 
of storm decay by exhibiting a warming in the vicinity of the over- 
shooting tops thresholds. Additional breakdown of the system in 
general is depicted by the remaining DCVMs (Figs. 22a-d) . The final 
GOES-1 Mg EIR image at 0400 GMT (Fig. 23, which should be compared 
with Fig. 19) displays the weakened, steady state of the activity in 
East Texas. 

At this point, it may prove meaningful to suggest that the DCVMs 
are also valuable when locating a collapse of an overshooting tops 
condition. As noted above, not only does the DCVM quickly establish 
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Fig. 22a. DCVM at 03/0300 GMT. 
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Fig, 22b. DCVM at 03/0315 GMT. 


52 


/-I 





Fig. 22c. DCVM at 03/0330 GMT. 
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the coldest region of the cloud, but it also clarifies any question 
that an observer might have when trying to determine whether he is 
looking at a light shade of gray (as in segment 5) corresponding to 
a much warmer temperature, or a colder portion (segment 8) of the Mg 
EIR curve. To illustrate this argument, compare Figs. 19, 21, and 23. 
Without similiar time frame DCVMs (Figs. 20a, 22a, and 22d) one could 
easily misinterpret the shades of gray and the conditions that are 
occurring in the atmosphere. The DCVMs prove that the count values 
corresponding to Fig. 23 are not as high (cold) as those of either 
Fig. 19 or Fig. 21. 

d. The Satellite-Derived Rainfall Estimations 

As mentioned earlier, the digital count value 208 was selected 
by empirical methods as the boundary threshold to indicate convective 
precipitation. The count values less than 208 were, therefore, 
considered as areas having little or no precipitation and were not 
used in this objective study. Since the Scofield and Oliver (1977) 

Mg EIR imagery technique uses temperature threshold values rather than 
digital count values as in this investigation, the application of 
their rainfall scheme required some adjustment. Naturally, this 
adjustment assumes that the overshooting tops criterion is expected 
or has already occurred. The basic change that was made was to expand 
Scofield and Oliver's "step 4" (estimates applying to the two coldest 
contours) of their rainfall estimation scheme. 

Since the rainfall estimation scheme using GOES imagery uses the 
two coldest contours in a particular thunderstorm system, the major 
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difficulty that a user must overcome is the identification of these 
coldest contours. The contours displayed may depend also on the 
character of the enhancement curve. This is where the digital counts 
are most readily usable because the highest (coldest) values are 
identified almost effortlessly, thereby enabling the user to pinpoint 
more quickly the area of maximum rainfall. 

Fig. 24 lists the rainfall amounts in mm %h ^ that were applied 
in using the cloud growth category versus digital count values in this 
study. Each count value in the left-hand column relates to a specific 
cloud top equivalent blackbody temperature. The next three columns 
relate the growth of the major axis of digital count contours between 
specific intervals of the DCVMs. These correspond to anvil growth 
rates that Scofield and Oliver used in their rainfall estimation 
scheme for %-h intervals. The rainfall amounts in mm %h ^ also were 
adapted from their imagery technique. The ^-h amounts were converted 
to mm ^h ^ for digital count values corresponding to their temperature 
thresholds. Appendix A outlines how the original breakdown of digital 
count values versus mm %h as displayed in Fig. 24, was devised. 

The satellite-derived rainfall estimates were arrived at by 
using an objective approach. The first decision made was to select 
the manner by which each hour's DCVMs would be handled. Again, by 
empirical methods, it was determined that the first four DCVMs of the 
individual hour would give an adequate representation of the hour's 
rainfall. This approach was used to keep the method as simple and 
straightforward as possible and also to prevent data overlap, which 
could possibly introduce some additional small errors. 
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Digital 

Count 

Value 

Increase 

in major 

axis of Anvil 

Digital Count 
Decreased to 
the 208 count 

Value Same or 
not less than 
value 

Decayed to less 
than 208 count 
value 

(along major axis 
>75 km 38-75 

of contour) 
km ^ 37 km 

208 

11.5 

7.5 

3.0 


2.0 


Trace 

209 

11.5 

7.5 

3.0 


2.0 




210 

11.5 

7.5 

3.0 


2.0 




211 

11.5 

7.5 

3.0 


2.0 




212 

11.5 

7.5 

3.0 


2.0 




213 

11.5 

7.5 

3.0 


2.0 




214 

13.5 

8.5 

4.0 


3.0 




215 

14.0 

8.5 

4.0 


3.0 




216 

14.5 

8.5 

5.0 


3.0 




217 

15.0 

8.5 

5.0 


4.0 




218 

15.5 

8.5 

5.0 


4.0 




219 

16.0 

9.0 

6.0 


4.0 




220 

16.5 

9.0 

6.0 


4.5 




221 

17.0 

9.0 

6.0 


4.5 




222 

17.5 

9.0 

6.0 


4.5 




223 

18.0 

9.0 

6.0 


5.0 




224 

18.5 

9.0 

6.0 


5.0 




225 

19.0 

9.0 

6.0 


5.0 



r 


Fig. 24. Rainfall in mm ' based on the growth or decay of the area of digital count value 
(adapted from Scofield and Oliver, 1977). 



Basically, the rainfall estimation scheme as developed from Fig. 


24 for the first hour is as follows: 


( 1 ) 


( 2 ) 


(3) 


(4) 


The 0100 GMT DCVM. 

(a) 3 mm %h ^ inside 208 contours (growth less than 37 km) 

(b) 4 mm %h ^ inside 214 contours (growth less than 37 km) 

The 0115 GMT DCVM. 

(a) 2 ram ^h ^ inside 208 contours (no change) . 

(b) 8.5 mm %h ^ inside 214-215 contours (growth 38-75 km). 

The 0130 GMT DCVM. 

(a) 2 mm %h ^ inside 208 contours. 

(b) 3 mm %h ^ inside 214-215 contours (no change) . 

(c) 5 mm %h ^ inside 217-218 contours (growth under 37 km) 

The 0145 GMT DCVM. 

-1 

(a) 2 mm %h inside 208 contours. 

(b) 3 mm %h ^ inside 214-215 contours (no change) . 

-1 

(c) 4 ram ^h inside 217-218 contours (no change) . 

(d) 6 mm %h ^ inside 220 contours (growth less than 37 km) 


This yields a total estimate for the hour of 0100-0200 GMT in mm. 


After the amount of precipitation was determined for each DCVM, 
the next step was to add them graphically. Before each analysis was 
concluded, it was necessary to displace the estimates an average of 
12 km toward the satellite subpoint to correct for the apparent 
location of the cloud versus the actual position of the cloud due to 
earth curvature. This procedure is recommended by Weiss (1978) and is 
outlined in Appendix B. Fig. 25 represents the graphic addition of 
DCVMs for 0100, 0115, 0130, and 0145 GMT by use of the estimation 
scheme outlined above. The same methods were used to determine Fig. 

26 and Fig. 27, the second and third hour, respectively. 
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Fig. 25. The satellite-derived rainfall estimation analysis for 
the first hour (0100-0200 GMT). 
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Fig, 26. The satellite-derived ra 
the second hour (0200-0300 GMT) . 
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Fig. 27. The satellite-derived rainfall estimation analysis for 
the third hour (0300-0400 GMT) . 
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The first and third hour rainfall estimations were quite simple 
and straightforward. The second hour estimation analysis, however, 
was somewhat complicated. The 0230 GICT DCVM depicts the case of 
extreme overshooting tops. This is a situation where an amplifica- 
tion factor of 6.5 mm %h ^ (as Scofield and Oliver endorse) must be 
incorporated into the analysis to account for the increased rainfall 
that is usually observed under such conditions. This situation is 
handled so that the top digital count value contour (in this case, 

222) is determined as the area that receives the 6.5 mm \h ^ of 
additional rainfall. This criterion was adapted to account for an 
anvil increase greater than 37 km when the digital count values also 
increase 5 or more counts (e.g., from 217 to 222) over a 30-min 
interval at levels at and above the tropopause. The 0230 GMT DCVM was 
the explosive period that fulfilled this criterion since the tropo- 
pause was approximately at the 212 digital count level and the over- 
shooting tops pushed 10 counts higher (10° C colder) as they increased 
5 counts during the 30-min interval prior to 0230 GMT. If the highest 
count level fulfilled all requirements except the 37-km anvil increase, 
the next lower count contour would be used if it was still above the 5 
count increase. 

The final step in this phase was to add graphically the indivi- 
dual hourly analyses into a 3-h total and produce a satellite-derived 
rainfall estimation analysis for this portion of the storm (Fig. 28). 
Table 4 depicts the amounts of rainfall derived at actual stations via 
this rainfall scheme for each hourly interval and the 3-h total. 
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Fig. 28. The 3-h total of satellite-derived rainfall estimation 
analyses (0100-0400 GMT) . 
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Table 4. Hourly rainfall estimates and 3-h total in mm 



Location 

First 

hour 

Second 

hour 

Third 

hour 

3-h total 

Belton Dam 

9 

9 

0 

18 

Conroe 

0 

2 

0 

2 

Eagle Lake 

0 

2 

2 

4 

Groesbeck 

2 

27 

10 

39 

Houston-WSMO 

0 

2 

0 

2 

Houston-Addicks 

0 

2 

0 

2 

Hou s t on- S at suma 

0 

2 

0 

2 

Jewett 

2 

10 

11 

23 

Lexington 

9 

29 

2 

40 

Love lady 

0 

6 

8 

14 

Post Oak School 

6 

11 

3 

20 

Red Rock 

10 

11 

0 

21 

Richmond 

0 

2 

0 

2 

Somerville Dam 

6 

10 

4 

20 

Stillhouse Hollow Dam 

9 

9 

0 

18 

Waco Dam 

9 

9 

0 

18 

Washington State Park 

0 

8 

4 

12 

Wheelock 

0 

24 

7 

31 
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e. Comparisons 

To compare properly the observed rainfall patterns to the 
satellite-derived rainfall patterns, it is important to relate the 
subjectivity of the raingage analyses. The raingage analyses were 
accomplished in a subjective manner inasmuch as the major part of the 
study area was without a well-defined, high-resolution, raingage 
network. The basic field, as outlined earlier, was analyzed from 18 
separate hourly recording stations. Even though these analyses were 
approached as objectively as possible, they will always be somewhat 
biased towards the individual performing the analysis. Furthermore, 
the 3-h total analysis (Fig. 13, p. 34) was performed similiarly to 
the hourly analyses with two additional restraints. The first 
restraint was that the 3-h analysis pattern had to fit graphically 
the addition of the individual analyses when they were combined, and 
the second, that the 3-h total rainfall pattern could not exceed the 
24-h analysis pattern (also derived subjectively, but with a greater 
number of reporting stations). On the other hand, the satellite- 
derived estimations are as objective as possible, within a range of 
certain reasonable restraints (Fig. 24, p. 58) to insure a field 
representative of the observed. 

With these facts in mind, comparisons can now be made of the 
hourly rainfall segments and the 3-h total portion of the storm. To 
maintain unbiased data, the correlations are based on the hourly 
rainfall station data versus the satellite-derived estimations at 
each recording station rather than a comparison of the subjective hand 
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analyses. The individual hourly raingage analyses and their respec- 
tive satellite-derived rainfall estimates will be discussed first, 
followed by the comparisons between the 3-h total observed versus 
estimated data. 

The first hour's data (Fig. 10, p. 31 and Fig. 25, p, 60) show 
similiarity with respect to the orientation of the rainfall pattern. 

Of particular interest are the first and third hour's rainfree areas 
and how they nearly match the observed when compared to the estimated. 
The observed versus estimated correlation for the first hour was 0.57 
for the 18 data pairs of information taken from Tables 2 and 4 (p. 29 
and p. 65). This was as expected because the rainfall amounts were 
low, which tended to allow the rainfall pattern to be more variable. 

The second hour received the most precipitation and the satellite 
estimates substantiate this (Fig. 11, p. 32 and Fig. 26, p. 61). Since 
the main thrust of this investigation is toward flash flood-type 
situations, it was important to discover the location of heavy precip- 
itation areas. The satellite estimation method exhibits ability for 
accomplishing this task. Two distinct regions of precipitation 
totaling 30 mm or greater are seen in both the observed and estimated 
analyses. Both analyses are in good agreement with their general 
rainfall patterns, and their 18-station data pair correlation was 0.85. 
This supports the statement that was alluded to in the above paragraph, 
that heavier rainfall will be more reliably indicated in the estima- 
tion technique, which will potentially provide more usable information 
to operational forecasters. 

The third hour (Fig. 12, p. 33 and Fig. 27, p. 62) was much like 
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the first, with less precipitation and a fairly widespread rainfall 
pattern. The northeast quadrant of the estimate depicts rainfall 
similiar to that observed. In the Post Oak School and Somerville Dam 
areas, however, it misses by over 10 mm. The third hour correlation 
was 0.56, about as expected. 

The 3-h total rainfall data sets will now be discussed. These 
data estimates should be considered separately from the hourly com- 
parisons and correlations because they were developed from those sets 
and remain dependent upon them. If the 3-h interval (0100-0400 GMT) 
was considered as the amount of time the "storm" took to cross the 
study area, then all the rainfall that was estimated would represent 
the storm's total. It is in this perspective that the 3-h total 
should be regarded. The data that covered the life of the storm would 
have the best rainfall patterns and correlations since the estimation 
technique is primarily based on four of the five DCVMs available for 
each hour. The small data loss, which was anticipated during each 
hour, will be recovered when the period of the storm is considered and 
the rainfall estimates for each hour are added up. 

Consequently, a comparison of the 3-h observed and estimated 
analyses (Fig. 13, p. 34 and Fig, 28, p. 64) indicates that the rain- 
fall patterns are strikingly similiar in areal extent for amounts up to 
30 mm and somewhat more widespread for greater amounts. Even though 
it seems obvious, it is worthwhile to note that isolated point values, 
represented by the raingages, would not (in a general case such as 
this) be expected to report the maximum amount of rain from a convec- 
tive system. With this in mind we can continue our comparisons. The 
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raingage analysis places the 40-mm isohyets over a smaller areal 
expanse than the satellite estimates. Rainfall appears to be over- 
estimated between Groesbeck and Lexington. This is confirmed if 
Fig. 9 (p. 28), the 24-h rainfall, is used as the upper limit of the 
precipitation. A visual comparison of the raingage analyses for 24-h 
and 3-h, and the 3-h satellite-derived estimate shows that the 3-h 
raingage analysis could be expanded, and then the comparisons between 
the observed and estimated would improve somewhat. The 18- station 
data pair correlation was 0,90. An examination in and around the 
Houston area, located in the southeast quadrant, indicated that little 
or no convective- type precipitation fell before 0400 GMT. Therefore, 
if the four reporting stations (Houston WSMO, Addicks, Satsuma, and 
Richmond) located in the southeast portion of that quadrant 
were removed from the station data pairs the correlation would improve 
to 0.92. 

Even though the 3-h 40-mm estimated contours appear larger than 
the 3-h observed 40-mm contours, they do provide valuable information 
and help to pinpoint heavy rainfall areas with reasonable accuracy. 
Table 5 was developed from Tables 2 and 4 and is offered as a quick 
comparison between the observed (raingage) and estimated (satellite) 
rainfall and the correlations computed each time interval. 

As mentioned earlier in the supplemental data section, the TAMU 
radar total VIL maps for selected areas were compared visually to 
DCVMs having a similiar time frame. Fig. 29 shows the 0223 GMT total 
VIL map superimposed on the displaced 0230 GMT satellite DCVM. 
Henderson (1979) examined this same time frame and computed a 
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Table 5. Hourly and 3-h observed and estimated rainfall in mm. 



Location 

First 
Obs Est 

Second 
Obs Est 

Third 
Obs Est 

3-h Total 
Obs Est 


Belton Dam 

4 

9 

7 

9 

0 

0 

11 

18 

Conroe 

0 

0 

0 

2 

0 

0 

0 

2 

Eagle Lake 

0 

0 

0 

2 

5 

2 

5 

4 

Groesbeck 

3 

2 

36 

27 

5 

10 

44 

39 

Houston-WSMO 

0 

0 

3 

2 

0 

0 

3 

2 

Houston- Addicks 

1 

0 

1 

2 

13 

0 

15 

2 

Hou s ton- Sat suma 

0 

0 

1 

2 

0 

0 

1 

2 

Jewett 

0 

2 

8 

10 

16 

11 

24 

23 

Lexington 

3 

9 

33 

29 

5 

2 

41 

40 

Love lady 

0 

0 

0 

6 

14 

8 

14 

14 

Post Oak School 

2 

6 

1 

11 

17 

3 

20 

20 

Red Rock 

0 

10 

20 

11 

0 

0 

20 

21 

Richmond 

0 

0 

1 

2 

4 

0 

5 

2 

Somerville Dam 

3 

6 

0 

10 

13 

4 

16 

20 

Stillhouse Hollow Dam 

4 

9 

4 

9 

0 

0 

8 

18 

Waco Dam 

12 

9 

4 

9 

1 

0 

17 

18 

Washington State Park 

0 

0 

3 

8 

8 

4 

11 

12 

Wheelock 

4 

0 

11 

24 

3 

7 

17 

31 

18-station correlation 

0. 

57 

0. 

85 

0. 

56 

0. 

90 
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Waco • 
Dam 



Fig. 29. The 0223 GMT total VIL map superimposed on the dis- 
placed 0230 GMT DCVM. The radar VIL isopleths begin with the analyzed 
values of 10 to 19 for the first blackened contour, and thereafter 
alternate with white, in steps of -10, until the highest VIL values of 
50-plus were detected approximately 20 km west of Jewett. 
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correlation of 0.39 for the northwest quadrant and a correlation of 
0.40 for the southwest quadrant. The figure clearly shows that the 
fields do not match up perfectly, but they do demonstrate some sig- 
nificant relationship. Additional total VIL maps were not available 
to pursue these comparisons further at this time. McAnelly (1979), 
however, is currently conducting research with these goals in mind. 
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6 . 


CONCLUSIONS 


This study demonstrated that an objective satellite-derived 
rainfall estimation technique is feasible with use of digital GOES 
IR data. Although these results are based only on one case, it 
appears that this approach can pinpoint areas of heavy, as well as 
lighter, precipitation with fairly accurate success. The ability to 
determine precipitation patterns and amounts in data-sparse areas 
would be an asset to many operational meteorologists, besides provid- 
ing valuable information to support investigative, research programs 
regarding the water budget. 

One of the shortcomings of the Scofield and Oliver (1977) rain- 
fall estimation scheme is that it requires an experienced individual 
to interpret properly the GOES VIS and EIR imagery. This investiga- 
tion requires that certain restraints be met and then uses the digital 
GOES IR data, which have the distinct advantage of 24-h surveillance 
of weather systems, to develop the precipitation estimates. The 
Scofield and Oliver (1977) and Scofield (1978a, 1978b) methodology 
requires careful interpretation of amplification factors. By contrast, 
I believe that the high digital count value variations in time and 
space will reveal objectively those situations where intense precipi- 
tation is occurring (e.g,, merging cell and intersecting thunderstorm 
lines). Furthermore, since the overshooting tops condition is in and 
above the anvil top, it is conceivable that the upper- level wind 
effects for short time intervals, as in this study, can be somewhat 
neglected to obtain operational estimates of heavy rainfall. 
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One of the foremost obstacles that must be surmounted is the 


reduction of the vast quantities of data to suitable spatial resolu- 
tion, This examination used McIDAS data, which have implied 4 km by 
4 km resolution, to demonstrate relatively high resolution rainfall 
patterns with areal variations that closely resemble the observed 
analyses. This tends to strengthen the argument even more for this 
kind of approach. 

The NWS currently is using the Scofield and Oliver method for 
estimating rainfall from hourly GOES Mg EIR imagery. Not only do they 
sacrifice resolution because of the mode of transmission of the 
imagery and the need to conserve manhours, they also may miss precip- 
itation amplifiers that could very likely be sandwiched between the 
hourly intervals of the EIR data. The technique used in this study, 
on the other hand, approaches the state of the art without forfeiting 
additional manhours. At the same time it will vastly improve and 
resolve the widespread precipitation pattern that is presently based 
on hourly GOES Mg EIR imagery. 

The DCVMs that were visually compared with the total VIL maps 
that were available indicated some fair correlations. One of the 
major problems could possibly be with the gridding, but I believe 
that a better understanding of cloud processes is necessary to reveal 
the true relationships between a VIL map and a DCVM. 

Although this single case does not prove conclusively that the 
methodology used herein is without flaws, it does demonstrate 
reasonable accuracy with a correlation of 0.92 for the convective 
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study period of the storm. This type of correlation accounts for 
approximately 85% of the variability in the data. This leads to a 
final conclusion that an accurate objective rainfall estimation 
technique based on digital geostationary satellite data is feasible. 
Furthermore, it could be one that should be developed further for 
operational use by meteorologists and hydrologists responsible for 
flash flood prediction. It could spell the difference between 
adequate warning and disaster. 
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7. 


RECOMMENDATIONS 


The following recommendations for use of digital GOES IR data 
for further satellite-derived rainfall estimation studies and flash 
flood detection investigation are presented: 

(1) This rainfall estimation approach should be tested against 
we 11“ documented rainstorms (e.g., flash floods such as those 
at Big Thompson Canyon, Johnstown, and Kansas City) to 
enable formulation of more conclusive results. These 
studies should also attempt to account for the rainfall 
underestimations made by Scofield (1978b) for the point 
maximum of the r a ingage measurements. 

(2) Additional research efforts should incorporate digital VIS 
data which have 1 km by 1 km resolution. With the use of 
the McIDAS gridding program devised by Henderson (1979), 
even higher resolution rainfall patterns may be possible. 
Furthermore, the comparisons made from the two types of 
digital data (VIS and IR) could reveal interesting 
conclusions . 

(3) Due to the areal limitations placed on the McIDAS data used 
in this study (240 km by 240 km), Atmospheric and Oceano- 
graphic Information Processing System (AOIPS) data should 
be utilized to conduct sirailiar research. The main advan- 
tage of the AOIPS is that it covers vast expanses of the 
United States, thus allowing studies in different regions 
from the same data base. 
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( 4 ) 


Consideration should be given to expanding segment 8 (the 
repeat gray levels) of the Mg EIR curve and then enlarging 
suspect flash flood areas on the Mg enhanced GOES imagery. 
This could be applied to the Scofield and Oliver (1977) 
rainfall scheme for pinpointing heavy precipitation areas 
while digital rainfall estimation schemes are being 
developed. 

(5) An algorithm should be developed that will allow adjust- 
ment, if necessary, to different geographic regions that 
are under surveillance. Development of the scheme should 
be incorporated into a computer program as a first step 
toward operational use. 

(6) Warm and cold (low and high digital count variations) anvil 
areas should be investigated to enable comparisons that 
would lead to a better understanding of the variations of 
the rainfall patterns associated with them. 

(7) Total VIL maps and DCVMs should be rigorously compared to 
examine their relationship and to enable understanding of 
the cloud processes taking place. 

(8) A mini-computer program incorporating the topography of a 
hydrological region (e.g., water shed) should be utilized 
with the rainfall estimates derived from various means 
(subjectively and/or objectively) to aid the meteorologist 
and/or hydrologist in operational flash flood prediction. 
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APPENDIX A 


This appendix outlines the approach used to relate the digital 
count values to the rainfall rate in mm %h Most of the information 
herein was adapted for use with digital GOES-1 infrared data and is 
based on the rainfall estimation scheme of Scofield and Oliver (1977)* 

When the attempt was made to adjust step 4 of the Scofield and 
Oliver decision tree, it became obvious that various variables would 
have to be converted to make them usable for this study. These 
included the rainfall rates versus various conditions of cloud growth, 
the Mg enhancement curve segments, and finally the cloud growth rates. 
Each of these will be discussed and then the adaptation of step 4 of 
the Scofield and Oliver decision tree will be presented. 

The rainfall rates were converted to mm %h ^ from in %h This 

was a simple task that required careful conversions into mm for each 
of the segments and various cloud growth conditions. After some 
experimentation with the digital data had taken place, it was 
determined that some small adjustments downward were necessary for 
the ”no change" (or count value not less than 208) category. Basi- 
cally, however, very little compensation such as this was made to the 
rest of the breakdown. 

Next, the Mg enhancement curve segments, which are given as 
temperatures only, were converted to digital count values. Since 
each segment of the Mg curve covered more than one digital count value, 
it was necessary to determine just how to relate the count values to 
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rainfall rates at various cloud growth stages. It was decided to 
spread the digital count values over that particular segment of a 
specified rainfall rate endorsed by Scofield and ^Oliver in step 4 of 
their estimation scheme. For instance, if the segment had the 
equivalent of 4 digital counts (e.g., varied over 4®C) and the rain- 
fall rate (for a particular cloud growth rate) increased by 2 mm %h ^ 

from the last segment, each digital count rainfall would be equivalent 

-1 “1 
to 0.5 mm %h . By the same token, if the change was 4 mm %h , each 

digital count value at that cloud growth rate level would represent 

1 mm %h”^. The digital count value versus rainfall rate breakdown 

was devised following this logical procedure. 

The final variable of interest was the cloud growth rate and 
how it related to the digital count values. Since the imagery scheme 
developed by Scofield and Oliver produced fairly good results, these 
same limits for cloud growth (anvil growth) were maintained for this 
investigation (e.g., an expansion in size of 2/3°Lat, approximately 
75 km in % h) . Some adjustment of this variable may be required in 
further studies. 

Fig. 30 lists the breakdown of digital count values versus rain- 
fall in ram %h ^ and their variance with respect to cloud (anvil) 
growth rates. These variables were adjusted as little as possible and 
were adapted almost directly from the imagery scheme of Scofield and 
Oliver. A quick comparison of Fig. 30 with Fig. 24 (p. 58) will 
disclose that the rainfall rates for the digital count values from 
208-213 were averaged for the cloud growth columns. This was done to 
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f 

Digital 

Count 

Value 

Increase in 
(along maj 
>75 km 

maior axis of Anvil 
or axis of contour) 
38-75 km ^ 37 km 

Digital Count Value Same or. 
Decreased to not less than 
the 208 count value 

Decayed to less 
than 208 count 
value 

177-186 

3.0 

2.0 1.0 

0.5 

Trace 

187-197 

6.0 

2.5 2.0 

0.5 



198-199 

7.0 

3.0 2.0 

0.5 



1 200-201 

8.0 

4.0 2.0 

1.0 



202-203 

9.0 

5.0 2.5 

1.0 



204-205 

9.5 

6.0 2.5 

1.0 



206-207 

10.0 

7.0 3.0 

1.0 



208 

10.5 

7.5 3.5 

2.0 



209 

11.0 

7.5 3.5 

2.0 



210 

11.5 

7.5 3.5 

2.0 



211 

12.0 

8.0 3.5 

2.0 



212 

12.5 

8.0 3.5 

2.0 



213 

13.0 

8.0 3.5 

2.0 



214 

13.5 

8.5 4.0 

3.0 



215 

14.0 

8.5 4.0 

3.0 



216 

14.5 

8.5 5.0 

3.0 



217 

15.0 

8.5 5.0 

4.0 



218 

15.5 

8.5 5.0 

4.0 



219 

16.0 

9.0 6.0 

4.0 



220 

16.5 

9.0 6.0 

4.5 



221 

17.0 

9.0 6.0 

4.5 



222 

17.5 

9.0 6.0 

4.5 



223 

18.0 

9.0 6.0 

4.5 



224 

18.5 

9.0 6.0 

4.5 



225 

19.0 

9.0 6.0 

4.5 




Fig. 30. Rainfall rates in mm based on the growth or decay of the area of digital count value 

(adapted from Scofield and Oliver, 1977). This figure was the basis for Fig. 24 (p. 58), 


allow an easier analysis of the DCVMs. When computer techniques are 
developed for this rainfall estimation scheme using digital data, 
even more expanded versions of Fig. 30 will be possible. 

Readers with further interest or questions regarding the 
technique used here may wish to refer to the paper by Scofield and 
Oliver (1977) for a more concise explanation of their work. 
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APPENDIX B 


CLOUD TOP DISPLACEMENT CORRECTIONS TCWARD SATELLITE SUBPOINT 

The appendix is provided for readers not familiar with the 
necessity to correct cloud top positions on satellite imagery. All 
information present below is from Weiss (1978) . 

Accurate positioning of cloud features on GOES imagery is 
dependent upon the height of the cloud tops and their distance from 
the satellite subpoint (SSP) . The apparent cloud location and the 
actual location coincide only over the SSP. Otherwise, the distance 
from the SSP increases due to earth's curvature and nadir angle. Fig. 
31 illustrates this problem. Therefore, for an individual to locate 
accurately a cloud on GOES imagery, he must apply a correction for 
earth curvature and nadir angle to obtain the appropriate displace- 
ment toward the SSP. 

The dimensions used for determining the cloud displacements are 
shown in Fig. 32, Angle 5^g is the great circle arc from the SSP to 
the cloud's apparent location on the image, h is the height of the 
cloud above the surface of earth, R is the mean earth radius. Re is 
the mean equatorial radius, A is the mean GOES altitude, and d is the 
location correction (the distance measured along the curved surface 
of earth). The additional angles in Fig. 32 ( ^ , cT » ^ » and ^ ) are 
utilized to calculate d. 
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Fig. 32. Dimensions used in calculating cloud location 
corrections (after Weiss, 1978). 


86 




To clear up possible confusion dealing with the mathematics, the 
sample calculation for a great circle arc of approximately 37® (the 
approximate interval from the SSP to CLL) is offered. 

R = 6371 km 
Re = 6378 km 
is = 37 ° 

A = 35, 793 km 
h = 13 km 

For triangle we can write 

R - (Re 

tan %(/^ _ 2T ) = ^ — — tan % (/^ + 

' R + (Re + A) 


= -0.737506 tan % (1) 

Since the sum of interior angles of a triangle equals 180°, we 
have for triangle 

/^ + er+ = 180° (2) 

and for + 3^ , we have ^ 

= 180° - = 143° (3) 

By substituting (2) and (3) into (1), we can write 

tan %(/^ - 'l^) = -0.737506 tan ^(180° - <(>^). (4) 

We are given that = 37°, so 

tan %(/5 - ) = -2.204172. (5) 

Evaluating (5) gives us 

^ - jr = -131.193826° . (6) 

Now, from (3) we have 
-vX = 143° . 
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By adding (3) to (6) we get 

/ _ ar = -131.193826“ 

+ ^ = 143.000000° 

11.806174“. 

Therefore 

/= 5.903087“. (7) 

By solving (3) for , we have 

er= 180“ ~ ((5g +/S ) - 137.096913“. (8) 

The law of sines for triangle gives 

h + R ^ R . (9) 

sin sin S 

Solving (9) for S yields 

sin S = 0.67937^ 

S = arcsin (0.679374) = 42.794750“. (10) 

For triangle , sxjnnning the interior angles we have 

S+ ^ = 180° . (11) 

Solving (ll)ifor-y gives 

nr^ = 0.108337“. (12) 

Finally arc length d is given by 
^ R'J? 

d = = 12.046471 km. 

Thus, for average heights of 13 km at a great circle arc of 
approximately 37“ , the displacement correction of the cloud tops is 
12 km toward the SSP. 
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